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The micro/nano encapsulation technology has acquired considerable attention in the fields of 
drug delivery, biomaterial engineering, and materials science. Based on recent advances in chemical 
particle synthesis, we propose a primitive model of an encapsulation system produced by the self- 
assembly of Janus oblate spheroids, particles with oblate spheroidal bodies and two hemi-surfaces 
coded with dissimilar chemical properties. Using Monte Carlo simulation, we investigate the en- 
capsulation system with spherical particles as encapsulated guests, for different densities. We study 
the anisotropic effect due to the encapsulating agent's geometric shape and chemical composition 
on the encapsulation morphology and efficiency. Given the relatively high encapsulation efficiency 
we find from the simulations, we believe that this method of encapsulation has potential practical 
value. 



The subject of encapsulation has received considerable 
attention in various fields of current research. For in- 
stance, in drug delivery, encapsulation with polymer mi- 
celles has some remarkable pharmacokinetic advantages, 
such as reduced systemic toxicity, reduced patient mor- 
bidity, and increased bioavailability and absorption effi- 
ciency. Encapsulation systems can also be used in spe- 
cific drug targeting in cancer treatment and controlling 
the drug release time. This has extensive therapeutic use. 
[T]-[6] In biomaterial research, the encapsulation of pro- 
teins could provide the building blocks for the formation 
of biodegradable scaffolds in tissue engineering. jTHTO] 
Furthermore, a polymer-encapsulating system promises 
to be industrially important for self-healing material. 
[TTJ [12] Several studies have been devoted to address- 
ing the issue of efficient, uniform encapsulation. Self- 
directed assembly provides a reasonable method through 
a bottom-up approach in micro/nanoscale systems. [13- 
[17J In particular, amphiphilic groups have been used 
in encapsulation since most encapsulated guests are 
hydrophobic molecules, that provide a high efficiency, 
are stable and have the prospect of broad application. 
[18l[19] 

In this article, we present a study of a new encapsu- 
lation method by means of the self-assembly of Janus 
oblate spheroids. Janus particles are two-faced particles, 
whose chemical makeup differs between the two sides. 
Due to the recent success in synthesizing these particles, 
there has been a significant amount of experimental and 
numerical work reported in recent literatures. [ 2QH22] 
For instance, Janus spheres have been studied in much 
of the experimental work, focusing on its synthesis and 
self-assembly properties. [23ti26] In addition, some theo- 
retical and numerical studies have also been carried out 
to investigate the phase diagram and cluster morphology 
based on simple, primitive, anisotropic potentials. [27]- 
[29j These studies show an anomalous thermodynamic be- 
havior in the gas-liquid coexistence region, due to the 
existence of stable micelle and vesicle aggregates in the 
gas phase (colloidal poor phase). [27 The inner hollow 



space and the exterior shielding (hydrophilic) surface of 
the shell structure associated with the micelle or vesicle 
aggregate makes it a good candidate to act as an encap- 
sulating agent. Furthermore, in the case of Janus oblate 
spheroids, it has two types of anisotropy; one is due to the 
Janus character of chemical composition and the other is 
due to its geometric shape. This reduces the rotational 
symmetry and also leads to the hollow space inside the 
shell structure formed by Janus oblate spheroids being 
larger than in the case of Janus spheres. These effects 
make it likely that oblate spheroids are more effective 
encapsulating agents than Janus spheres. Thus in this 
paper we investigate the performance of encapsulation 
by Janus oblate spheroids. We believe that Janus oblate 
spheroids may lead to a high efficiency of encapsulation 
and serve as a platform for the implantation of functional 
groups. 



MODEL AND SIMULATION METHOD 

In our model, the encapsulating agents are oblate 
spheroidal Janus particles, defined as a hemi- spheroid 
patchy model of oblate spheroids (ellipsoids with the 
lengths of the principal axes denoted byc<a = 6 = cr/2 
and aspect ratio e = c/a); the encapsulated guests are 
chosen to be isotropic hydrophobic spheres with diameter 
R = a. There are three kinds of inter-particle interac- 
tions, ellipsoid-ellipsoid (/7ee), ellipsoid-sphere (Ues), and 
sphere-sphere (Uss)^ respectively. 

Analogous to the patchy spherical model introduced 
by Kern and Frenkel, we choose the ellipsoids to interact 
through a pair-potential that depends on their separation 
and orientation: Uij = U f{fij ^hi^hj). [30] The standard 
square-well potential has been used in earlier work for 
Janus spheres; however, the determination of the accu- 
rate spatial relation between ellipsoids is computationally 
time-consuming. Thus, as in |31j, we introduce a 'quasi- 
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where uq is the well depth, H{x) denotes the Heaviside 
step function, rij is the center-to-center distance between 
particles and fij is the direction of particle separation 
vector fij = fi — fj^ fii^hj are the patch orientations 
which are along the semi- minor axes, as shown in Fig. [T] 
We treat isotropic hydrophobic sphere as a unique Janus 
ellipsoid, whereas its e = 1.0, patch covers its whole sur- 
face with orientation fii = fji{nj = fij). In this setting, 
it meets the criterion that the three types of interactions 
{Uee^ Uesi and Uss) could be unified in the potential form 



is introduced as an approximation to character- 



ize the spatial relation from representing molecular inter- 
actions by a Gaussian model potential, such that there is 



no overlapping interaction if Vij > ai 



Therefore, in 



our case H{aij + 0.5a — rij) represents a 'quasi- square- 
well' potential with width 0.5cr. For the hydrophobic 
interaction between spheres, the potential is a standard 
square- well potential with well-depth denoted as Uss and 
a cut-off range Tc = 1.5a (since aij = a). To summarize, 
a) the interactions between two hydrophobic parts (hemi- 
surfaces of Janus ellipsoids, or entire spherical surfaces) 
are 'quasi-square-well' potentials; b) the interaction be- 
tween a hydrophilic region and other regions is just a 
hard-core repulsion. 




FIG. 1. (color online). Illustration of Janus oblate spheroids 
(e = 0.4). The hydrophobic side is coded in red and the hy- 
drophilic side in blue. Red hemi- surfaces attract each other; 
blue hemi- surfaces interact with other surfaces through a sim- 
ple hard-core repulsion. 

We use standard Monte Carlo (MC) simulation in the 
NVT ensemble to study the encapsulation process of this 



system. We choose small three-dimensional (3D) systems 
of size L = 2Aa (as the simulation becomes quite com- 
putationally expensive with increasing particle number) . 
Periodic boundary conditions are enforced to minimize 
wall effects. All simulations start from a random initial 
monomer conformation and the results are averaged over 
tens of runs. In the following, a is taken as unit of length 
and uq the unit of energy. The temperature T is also ex- 
pressed in the unit of Uq (with the Boltzmann constant 
kB — 1). The other parameters are denoted as follows: 
the number of ellipsoids A/'g, the number of spheres A/'^, 
the ratio of particle numbers a = N^/Ns^ and the total 
number density p = {N^ -\- Ns)/L^. 



RESULTS AND DISCUSSION 

Encapsulation is defined such that each guest (sphere) 
is sufficiently surrounded by agents (oblate spheroids). 
To be specific, in our model, the self-assembly of the 
Janus ellipsoids leads to aggregates with one or more 
outer shells; the outmost surfaces are (mostly) hy- 
drophilic. These shells surround hydrophobic centers 
consisting of one or more spheres. To study the mech- 
anism of the encapsulation as well as the efficiency de- 
pendence on the ellipsoid aspect ratio, we need to find 
a way to quantify the degree of encapsulation that char- 
acterizes how well each sphere is surrounded by ellip- 
soids. In the first case we study, we choose a low density, 
p = 0.02 and temperature T = 0.25, since we find that 
this thermodynamic point is in the gas phase (micelle) 
region. We choose a ratio of particle numbers a = 10, 
as this equals the number of particles in a stable micelle 
of Janus spheres [27 and set Uss = so that there is no 
hydrophobic attraction between spheres (just hard-core 
repulsion). We choose Ugs = initially to make it fa- 
vorable that only one sphere is encapsulated. Later we 
consider the case for which Ugg ^ 0. We carry out simu- 
lations up to 8 million Monte Carlo steps (MCS). From 
the results (not shown here) , we find a window of simula- 
tion time around 2.5 million MCS when mono-dispersed 
clusters form (with mostly single guests inside), which 
co-exist with a few monomers and small oligomers. All 
the simulation results presented in this paper are from 
the runs in this window, at 2.5 million MCS. 

We investigate the clusters formed by plotting the dis- 
tribution of the particle distance from the cluster center 
of mass. A typical example is shown in Fig. [2] From 
the result, we find that inside this cluster there are 26 
particles with one sphere (colored green) and mostly two 
layers (shells) of ellipsoids; the inner layer consists of 12 
ellipsoids. The number 12 is typical of what we find for 
the (inner) layer of ellipsoids in other clusters and is also 
reminiscent of the size of stable Janus sphere micelles. 
[29 Based on this observation, we define the 'ideal en- 
capsulation' as corresponding to the case in which each 
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TABLE I. Efficiency (%) of 'ideal encapsulation' at tempera- 
ture T = 0.25 



^ \ 


1.0 0.7 0.6 0.5 0.4 0.3 


0.01 
0.02 
0.04 
0.06 


15.5 9.1 13.6 10.9 7.3 11.8 
24.7 21.3 27.3 28.0 24.0 21.3 
28.9 33.3 44.4 35.0 30.0 26.7 
23.3 26.7 25.9 28.9 32.9 40.7 



sphere is surrounded by at least 12 ellipsoids so that it 
forms a cluster with a coniplete shell. Then we calculate 
the efficiency 77 by dividing the number of spheres that 
are 'ideally' encapsulated (A^|) by the number of total 
spheres inside the system (A/'g), i.e., r] = N^/Ns. 
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FIG. 2. (color online). Plot of distance r{i) of the i^^ particles 
to the center of mass of the cluster, corresponding to the 
configuration shown in the inset. 

We choose a = 12,T = 0.25, and run the simulations 
for different densities and aspect ratios. We show the 
results in Table [l[ The table shows that: a) as the aspect 
ratio e varies from 1.0 to 0.3, the efficiency increases and 
achieves its peak value at e = 0.6, then it decreases; b) 
the same trend is found in the variation with density and 
reaches another peak at p = 0.04. Two typical runs for 
density p = 0.01 and 0.06, with an ellipsoid aspect ratio 
e = 0.4, are presented in Fig. [3j 

The behavior shown in Table [J can be understood 
in terms of the balance between enthalpy- driven and 
entropy- driven effects. To illustrate this competition, 
we first examine the numerical results for the second 
virial coefficient, that characterizes the pair inter- 
action between the particles, for different aspect ra- 
tios at three temperatures. We plot B2 for the Janus 
ellipsoid-ellipsoid interaction, {B2)eei and a scaled B2 
of the Janus ellipsoid- sphere interaction, with (52)es = 
(52)es/(^2)ee7 as a fuuctiou of € in Fig. [4] The stars in 
Fig. [4] A represent the theoretical prediction for Janus 




FIG. 3. (color online). Morphology of the entire system for 
an ellipsoid aspect ratio e = 0.4. (A) Initial configuration of 
number density p = 0.01. (B) Encapsulation system (p = 
0.01) after 2.5 million MCS, efficiency ry = 9.1%. (C) Initial 
configuration of number density p = 0.06. (D) Encapsulation 
system {p = 0.06) after 2.5 million MCS, efficiency r] = 48.1%. 

spheres. [30l |3l] Fig. [4] shows a monotonic behavior of 
these quantities: a) {B2)eei the magnitude of the effective 
interactive strength increases as temperature decreases 
and the aspect ratio goes up. This means that particles 
combine more strongly at lower temperature; in addition, 
the larger the aspect ratio, the stronger is the attraction 
to each other (noticing that with same interaction range, 
interaction space is proportional to e); b) (^Des indicates 
how larger the Janus ellipsoid-sphere interaction is com- 
pared to the Janus ellipsoid-ellipsoid interaction. As one 
would expect, the greater this value is, the more likely an 
ellipsoid chooses to combine with a sphere instead of an- 
other ellipsoid. We see the advantage of a smaller aspect 
ratio ellipsoid from the curves. This also suggests that 
such an advantage becomes more important at higher 
temperatures. However, the absolute values of B2 are 
relatively small at high temperatures, i.e., the interac- 
tions are relatively weak, which is not helpful in forming 
stable encapsulation structures. From this discussion, we 
conclude that enthalpy-driven effect should lead to the 
efficiency increasing as the aspect ratio decreases. 

Next we evaluate the distribution of the scalar prod- 
uct ni • n2 for all pairs of bonded ellipsoid particles in 
three cases, e = 0.3,0.6 and 1.0 respectively, and give 
the results in Fig. |5]A. We illustrate the relative shapes 
and orientations of the bonded particles by the confor- 
mations shown there. For e = 1.0 (Janus spheres), the 
conformation of pair particles mainly lies in two regimes, 
facing each other and nearly parallel to each other. This 
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FIG. 4. (color online). Plot of second virial coefficient {B2) vs. ellipsoid aspect ratio at three temperatures, T = 0.25, 0.28, and 
0.33. (A) B2 of Janus ellipsoid-ellipsoid interaction, {B2)ee- Color stars indicate the values from theoretical prediction. (B) 
Scaled B2 of Janus ellipsoid-sphere interaction on the (B2)ee, i-e. (5J)es = {B2)es / {B2) ee- 



ls favorable in forming large clusters with complex struc- 
tures, providing two types of dissimilar building blocks. 
For e = 0.6 and 0.3, the curves are localized around single 
peaks. However, the orientation configuration for e = 0.3 
is more spread out, which indicates a larger entropy ef- 
fect due to the freedom of the arrangement of the patch 
direction. Also, the peak position of e = 0.3 around 
0.20 implies that the angle between two patch orienta- 
tions is Op = co5~^ (0.20) = 78.5°. Therefore, the angle 
formed between the equatorial planes of the pair of el- 
lipsoids Oe = 180° - 78.5° = 101.5°. Considering the ge- 
ometric configuration, the shell structure formed by the 
pairs with such an internal angle and the other analo- 
gous pieces might have a smaller hollow space inside. As 
a consequence, this is not favored in encapsulation. On 
the contrary, for e = 0.6, the curve is more localized and 
the peak dominates the distribution, with its pair config- 
uration open enough {0^ = 180° - cos-^{0.58) = 125.5°) 
that the structure aids in the self-assembly of encapsu- 
lation. This argument finds additional support in Fig. 
[5] B, which shows the cluster size {Np) distribution for 
three aspect ratios, e = 1.0 has a large and extended 
distribution (up to Np = 96), while the histograms of 
e = 0.3 and 0.6 are much compressed, e = 0.3 is peaked 
sharply around a small size, since it is easier to form 
a shell structure with a shielding outer surface within 
small clusters, considering that the pair configurations of 
e = 0.3 have smaller internal angles. Taking the two- 
dimensional case for example, the equiangular polygon 
with an internal angle leads to the number of sides 
n = 360°/(180° — 0). Hence the smaller the internal an- 
gle, the fewer sides and fewer vertices the polygon could 
have. In conclusion, the entropy-driven effect becomes 
obvious for ellipsoids with small aspect ratios and causes 



the reduction of encapsulating efficiency. As a conse- 
quence, it is reasonable to have e = 0.6 as a practical 
choice for Janus oblate spheroids in encapsulating guests 
through the self-assembly. 

Now we move on to discuss another factor that might 
affect the efficiency, i.e., the interaction between encapsu- 
lated guests {uss)- We start the new runs at a = 12, T = 
0.25 for different densities {p = 0.02 and 0.04) and aspect 
ratios (e = 0.5 and 0.6) for three values of Uss- We con- 
sider these more realistic cases in order to compare with 
practical applications, to shed light on potential experi- 
mental work. We introduce a new method to calculate 
the efficiency, by direct observation from different 3D per- 
spective views. For each run, we view all the clusters gen- 
erated individually via 3D rotation and count the total 
number of sphere (s), n^, not sufficiently surrounded to be 
considered encapsulated. We repeat these observations 
for each run three times, average the number of spheres 
'at large' n = (n^) and calculate the efficiency r]^ through 
{Ns — n)/Ns. The results are given in Table [H} From 
these values, we find relatively high efficiencies, which 
lends credence to the argument that the self-assembly 
of Janus ellipsoids provides a possible mechanism for ef- 
ficient encapsulation. Furthermore, we have not found 
any evidence that the interaction between encapsulated 
guests would change the encapsulation dramatically at 
low densities, for the cases we consider here. It is pos- 
sible that fluctuations and time variation in the cluster 
equilibration process of different sizes produce the dif- 
ferences in the efficiencies. However, this leaves an open 
question for the further study of this model at higher den- 
sities. To end this section, we show two typical snapshots 
of the final stages in Fig. |6]for e = 0.6. In both systems, 
only two guests (spheres) have not been entrapped sufii- 
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FIG. 5. (color online). Distribution of (A) orientation correlation over all pairs of bonded ellipsoid particles and (B) cluster 
size (Np) inside the system after 2.5 million MCS for different ellipsoid aspect ratios, e = 0.3,0.6, and 1.0 respectively. 



TABLE II. Efficiency (%) of encapsulating from observation 
at temperature T — 0.25 



e \ 


0.02 


0.04 


Uss/uq 




60.2 


62.5 





0.5 


73.6 


61.4 


0.5 




65.0 


65.9 


1.0 




80.7 


81.8 





0.6 


74.6 


63.9 


0.5 




75.9 


81.1 


1.0 



• 




ciently, i.e., only two are not encapsulated. In addition, 
we have also studied the cases in which there is no in- 
teraction between guests and agents, namely U^s = 0; in 
this case encapsulations were rarely found. The results 
suggest that by screening the interaction between guests 
and agents, one could switch off the encapsulation. In 
this case one would expect to see separate phase separa- 
tion for each of the two species. 



CONCLUSION 

In summary, we have proposed a potential method 
for efficient encapsulation of interacting guests, through 
the self-assembly of Janus oblate spheroids, which based 
on recent advances in chemical synthesis of Janus par- 
ticles. We have defined a model correspondingly and 
used Monte Carlo simulations to investigate the encap- 
sulation process. We find that under the competition 
of enthalpic and entropic effects, it is most favorable to 
use Janus oblate spheroids with an aspect ratio e = 0.6, 
at least in the range of ratio values varying from 0.1 to 
1.0 (via a step size 0.1). We also examine other possi- 
ble factors in determining the efficiency, such as density. 



FIG. 6. (color online). Morphology of entire system for ellip- 
soid aspect ratio e = 0.6 after 2.5 million MCS. (A) Encapsu- 
lation system p = 0.02, i^ss = uq. (B) Encapsulation system 
p = 0.04,1^35 0. 



temperature, and interaction strength etc. The results 
turn out to support the argument that a method based 
on anisotropic properties introduced by geometric shape 
and chemical composition could be quite general in its 
application. Thus we believe the methodology of encap- 
sulation we propose here could be of practical value and 
raise questions for future study. 
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